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H202  reduction  reaction  (HPRR)  is  studied  on  both  graphene  (GNF)  and  nitrogen  doped  graphene 
nanoflakes  in  0.1  M  Na2S04  solution  by  rotating  disk  electrode.  The  XPS  results  indicate  that  N-doped 
graphene  nanoflakes  with  high  nitrogen  content,  32  at%N  (N-GNF32),  are  synthesised  successfully  by  an 
inductively-coupled  thermal  plasma  (1CP)  reactor.  Pyridinic,  pyrrolic  and  graphitic  N  species  contribute 
up  to  67%  of  the  total  nitrogen.  Kinetic  parameters  such  as  Tafel  slope  and  stoichiometric  number  suggest 
that  HPRR  occurs  by  the  same  mechanism  on  both  GNF  and  N-GNF32.  Although  nitrogen  does  not  change 
the  mechanism  of  HPRR,  the  results  indicate  that  the  reaction  rate  of  H202  reduction  is  enhanced  on  N- 
GNF32.  The  exchange  current  density  of  H202  reduction  based  on  the  active  surface  area  of  N-GNF32  is 
(8.3  ±  0.3)  x  1CT9  A  crrr2,  which  is  6  times  higher  than  the  value  determined  for  GNF.  The  apparent 
number  of  electrons  involved  in  the  process  suggests  that  H202  decomposition  competes  with  H202 
reduction  on  both  catalysts.  Evaluation  of  the  apparent  heterogeneous  reaction  rate  constant  and  the 
Tafel  slope  indicate  that  simultaneous  reduction  of  02  and  H202  is  negligible  on  the  N-GNF32.  On  the 
other  hand,  the  reduction  of  02  and  H202  occurs  simultaneously  on  the  GNF  surface. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 
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Nitrogen  doped  carbon  nano-materials  have  been  studied 
widely  to  replace  the  conventional  metallic  catalysts  in  fuel  cells 
[1,2].  N-doping  is  a  method  to  increase  the  catalytic  activity  of 
carbon  nano-materials  toward  the  oxygen  reduction  reaction 
(ORR);  however,  the  phenomena  responsible  for  the  enhancement 


http://dx.doi.org/10.1016/j.jpowsour.2014.01.114 
0378-7753/©  2014  Elsevier  B.V.  All  rights  reserved. 


S.J.  Amirfakhri  et  al.  /  Journal  of  Power  Sources  257  (2014)  356-363 


357 


of  the  catalytic  activity  for  nitrogen-doped  carbon  nano-materials 
are  largely  unknown  and  still  controversial  [1,3,4], 

Nitrogen  doped  graphene  is  one  of  the  new  catalyst  structures 
that  has  received  a  great  of  attention  in  recent  years  particularly  in 
the  context  of  the  oxygen  reduction  reaction  (ORR)  [5,6],  Graphene 
is  a  two  dimensional  form  of  carbon  which  shows  special  charac¬ 
teristics  such  as  high  crystallinity  helping  to  maintain  stability,  high 
specific  surface  area  and  electrical  conductivity.  Different  types  of 
nitrogen  species  in  the  graphene  lattice  present  different  catalytic 
performance  toward  ORR.  For  instance  Lai  et  al.  [7]  showed  that 
while  the  graphitic  N  species  noticeably  increase  the  limiting  cur¬ 
rent  density,  the  pyridinic  N  shifts  the  onset  potential  towards  more 
positive  values.  Moreover,  pyridinic  N  species  might  change  the 
ORR  mechanism  from  the  series  (Eq.  (1))  to  the  direct  dominated 
pathway  (Eq.  (2)). 


02  +  2e~  +  2H+- 

-h2o2 

(1) 

02  +  4e~  +  4H+- 

■+2H20 

(2) 

Although  ORR  has  been  studied  extensively,  few  results  have 
been  reported  on  the  reduction  of  H2O2  using  N-doped  graphene  as 
a  catalyst.  Hydrogen  peroxide  has  been  investigated  as  an  alter¬ 
native  oxidant  to  oxygen  in  fuel  cells  for  several  reasons  [8—10], 
First,  H2O2  has  a  much  higher  solubility  than  oxygen  in  water;  it  is 
thus  easily  accessible  in  the  liquid  phase  at  higher  concentrations 
than  oxygen.  The  use  of  an  aqueous  oxidant  which  is  more  soluble 
in  water  than  dissolved  oxygen  is  a  possible  solution  for  an  increase 
of  the  power  density  of  fuel  cells.  Second,  the  hydrogen  peroxide 
reduction  reaction  (HPRR)  shows  a  lower  activation  barrier  in 
comparison  to  the  oxygen  reduction  reaction  (ORR).  The  HPRR  is  an 
electrochemical  reaction  involving  two  electrons  being  transferred, 
while  a  four-electron  transfer  reaction  occurs  in  the  ORR.  Finally, 
H2O2  is  a  suitable  candidate  for  space  and  underwater  applications 
where  air  is  not  readily  accessible  to  supply  oxygen  for  fuel  cells. 

The  most  important  disadvantage  of  using  H2O2,  however,  is  the 
H2O2  decomposition.  H2O2  can  be  reduced  following  Eq.  (3)  and/or 
decompose  to  oxygen  on  the  electrode  surface  through  Eq.  (4). 
H2O2  decomposition  causes  the  partial  loss  of  the  reactant,  de¬ 
creases  the  power  density  of  the  fuel  cell  and  can  also  block  the 
reactant  flow  inside  small  microfluidic  fuel  cells  [11,12], 

H202  +  2e+2H+^2H20  (3) 

2H202  — >02  +  2H20  (4) 

Among  the  noble  metal  catalysts,  gold  exhibits  high  catalytic 
activity  towards  H202  reduction  while  partially  inhibiting  H202 
decomposition  13,14].  Non-metal  as  well  as  non-noble  metal  cat¬ 
alysts  of  various  types  have  also  been  studied  for  HPRR  [15—18], 
Recently,  Shao  et  al.  [19]  studied  HPRR  on  N-doped  graphene 
qualitatively  in  a  quiescent  solution.  The  results  of  cyclic  voltam¬ 
metry  and  chronoamperometry  indicated  that  N-doped  graphene 
produces  much  higher  current  densities  in  comparison  to  pure 
graphene.  However,  no  information  has  been  reported  on  the  ki¬ 
netic  parameters  such  as  the  exchange  current  density  or  the 
mechanism  of  HPRR  occurring  on  N-doped  graphene  structures. 

Very  recently  the  H202  reduction  mechanism  has  been  studied 
by  Wu  et  al.  20]  on  both  graphene  and  N-doped  graphene  using 
density  functional  theory  (DFT)  calculations.  Their  results  indicate 
that  the  H202  adsorption  on  both  graphene  and  N-doped  graphene 
occurs  through  a  weak  physical  bond  (physisorption).  Moreover, 
two  paths  are  suggested  for  the  H202  reduction  on  both  graphene 
and  N-doped  graphene.  In  the  first  path  the  adsorbed  H202,  there¬ 
after  called  (H202)ads,  receives  an  electron  to  produce  adsorbed  OH, 


(OH)ads.  Then  (OH)ads  receives  another  electron  to  produce  water.  In 
the  second  path,  the  (H202)ads  first  breaks  to  two  (OH)adS  and  then 
(OH)ads  obtains  one  electron  to  produce  water.  In  the  latter  path  the 
electrochemical  step  is  repeated  twice  to  complete  the  process. 
According  to  the  calculation  of  the  Gibbs  free  energy  of  different 
reactions  involved  in  these  two  paths,  the  first  path  is  suggested  by 
Wu  et  al.  for  the  H202  reduction  on  both  graphene  and  N-doped 
graphene.  It  is  also  suggested  that  the  pyridinic  N  atoms  should 
exhibit  the  best  performance  for  the  H202  reduction. 

In  this  research,  the  catalytic  activity  of  nitrogen  doped  gra¬ 
phene  nano-flakes  (N-GNF)  with  32%  nitrogen  content  (thereafter 
labelled  N-GNF32)  and  pure  graphene  nano-flakes  (GNF)  are 
studied  experimentally  for  HPRR  in  a  neutral  electrolyte  (0.1  M 
Na2S04)  by  the  rotating  disk  electrode  (RDE)  technique.  Pyridinic, 
pyrrolic  and  graphitic  N  species  constitute  67%  of  the  total  nitrogen 
of  N-GNF32.  The  objectives  of  the  current  work  include  (1)  the 
determination  of  the  kinetic  parameters  for  HPRR  on  both  gra¬ 
phene  and  N-doped  graphene,  (2)  an  evaluation  of  H202  decom¬ 
position  according  to  the  determined  kinetic  parameters,  (3)  a 
comparison  of  the  reaction  rates  on  the  catalysts,  and  (4)  the  sug¬ 
gestion  of  a  mechanism  for  the  H202  reduction  on  both  graphene 
and  N-doped  graphene  according  to  the  kinetic  parameters  deter¬ 
mined  experimentally.  Finally,  the  produced  current  density  on 
GNF  and  N-GNF32  are  compared  with  gold. 


2.  Experimental 

2.2.  Materials  and  measurements 


The  synthesis  methods  of  GNF  and  N-GNF32  have  been  reported 
in  detail  elsewhere  [21].  In  summary,  an  inductively-coupled 
thermal  plasma  (ICP)  reactor  uses  the  dissociation  of  methane 
precursor  to  produce  pure  graphene  nanoflakes  through  controlled 
homogeneous  nucleation  of  the  carbon  vapour  at  high  tempera¬ 
tures.  N-doped  graphene  with  different  nitrogen  content  can  be 
generated  in  situ  using  a  nitrogen  plasma  as  a  post  treatment  on  the 
GNF  powders  within  the  reactor.  Transmission  electron  microscopy 
(TEM),  and  X-ray  photoelectron  spectroscopy  (XPS)  were  used  to 
characterize  the  morphology  and  evaluate  the  elemental  compo¬ 
sition  of  the  samples.  The  Brunauer— Emmett— Teller  (BET)  method 
from  nitrogen  gas  adsorption— desorption  isotherms  at  77  I<  was 
used  to  determine  the  specific  surface  area  of  the  powders.  The 
details  of  the  working  electrode  preparation  and  the  electro¬ 
chemical  tests  are  presented  in  the  Supporting  information,  section 
1  and  2.  The  H202  concentration  in  the  control  experiment  was 
measured  using  the  Idiometric  method  [22], 

The  active  surface  area  of  the  electrodes  was  determined  by 
voltammetry  method  [23],  In  this  method,  the  double  layer 
charging  current  (i'di)  is  measured  at  different  scan  rates.  A  linear 
relation  is  expected  between  the  double  layer  charging  current  and 
scan  rate  sr  =  d£/dt,  its  slope  being  used  to  determine  the  apparent 
total  double  layer  capacitance  of  the  electrode  (C/ \)  according  to  Eq. 
(5)  where  Q  is  the  charge  transferred. 


r  _  dQ  _  jdl  _  'dl 

A  d£  sr 


(5) 


The  active  surface  area  is  then  obtained  by  dividing  the  apparent 
total  double  layer  capacitance  to  the  double  layer  capacitance 
available  from  the  literature.  The  specific  active  surface  area  is 
determined  by  dividing  the  active  surface  area  with  the  catalyst 
mass.  The  exchange  current  densities  are  normalized  to  the  active 
surface  areas  in  order  to  determine  the  intrinsic  effect  of  the 
catalyst  structure  on  HPRR. 
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Table  1 


XPS  results  of  the 

samples(at%):  GNF  and 

N-GNF32. 

Sample 

N 

C 

O 

GNF 

0 

98.1  ±  0.3 

1.8  ±  0.4 

N-GNF32 

32.4  ±  0.2 

63.7  ±  0.5 

3.3  ±  0.5 

3.  Results  and  discussion 

HPRR  was  investigated  on  CNF,  N-GNF32  and  gold  in  a  neutral 
solution  (0.1  M  Na2S04)  by  RDE.  Before  the  tests,  the  stability  of 
H2O2  in  0.1  M  Na2S04  solution  was  evaluated  to  ensure  that  H2O2 
decomposition,  if  any,  only  occurs  on  the  electrode  and  not  in  the 
bulk  solution.  For  this  test,  a  concentration  of  0.1  M  of  H2O2  in  a 
0.1  M  Na2S04  solution  was  used  as  a  control  experiment,  the  H202 
concentration  being  measured  in  30  min  time  intervals.  The  results 
of  this  control  experiment  indicated  that  H2O2  concentration  star¬ 
ted  to  decrease  slightly  after  1  h.  After  2  h,  more  than  97%  of  H2O2 
was  still  available  in  the  solution  indicating  that  H2O2  is  reasonably 
stable  for  sufficiently  long  periods  of  time  in  this  solution  to  com¬ 
plete  the  electrochemical  tests. 

3.1.  Catalyst  characterization 

The  nitrogen,  carbon  and  oxygen  ratios  on  the  series  of  GNF  and 
N-GNF32  samples  used  in  the  present  experiments  are  given  in 
Table  1  following  XPS  analysis,  confirming  the  large  amount  of 
nitrogen  present  in  the  N-GNF32  [21].  The  results  demonstrated 
that  atomic  nitrogen  levels  over  30  at%N  on  the  carbon  support 


(a) 


were  achieved.  The  distribution  of  the  nitrogen  was  24%,  28%,  and 
15%  in  pyridinic,  pyrrolic  and  graphitic  sites,  respectively.  No  metal 
impurities  were  detected  in  the  synthesised  samples. 

The  TEM  images  in  Fig.  1  show  the  structure  of  both  GNF  and  N- 
GNF-32.  These  are  made  of  an  average  of  10  graphene  layers 
(number  of  layers  varying  between  5  and  20)  and  have  planar  di¬ 
mensions  typically  in  the  range  between  50  and  100  nm.  The  large 
nitrogen  content  is  believed  to  occur  mainly  on  the  surrounding 
edges  of  the  stacked  graphene  planes  structure. 

The  BET  specific  surface  area  of  the  samples  was  measured  as 
345  and  130  m2  g_1  for  GNF  and  N-GNF32,  respectively.  These 
values  are  much  lower  than  the  specific  surface  area  of  an  individual 
graphene  sheet  (2600  m2  g_1 )  suggesting  the  presence  of  multilayer 
graphene  and  N-doped  graphene  sheets  [24  .  This  result  is  consis¬ 
tent  with  the  observation  from  TEM.  Moreover,  the  BET  specific 
surface  area  of  the  N-GNF32  is  unexpectedly  3  times  lower  than  the 
GNF.  Generally,  the  specific  surface  area  of  graphene  is  expected  to 
increase  after  N-doping  because  of  inducing  more  disorder  in  the 
graphene  lattice.  Agglomeration  of  the  graphene  nanoflakes  in  the 
post-treatment  process  to  produce  N-GNF32  might  be  a  reason  for 
the  low  BET  results  determined  in  our  work  although  the  TEM  ob¬ 
servations  do  not  seem  to  show  any  morphological  changes.  A 
decrease  of  the  specific  surface  area  of  graphene  after  nitrogen 
functionalization  has  also  been  reported  by  Sheng  et  al.  25], 

3.1.1.  Determination  of  specific  active  surface  area 

Fig.  2(a)  shows  the  polarization  data  of  a  GNF  electrode  in  0.5  M 
H2SO4  solution.  Cyclic  voltammograms  in  the  potential  range  from 


Fig.  1.  TEM  images  of  the  catalysts:  GNF  with  low  (a)  and  high  magnifications  (b),  and  N-GNF32  with  low  (c)  and  high  magnifications  (d). 
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Fig.  2.  (a)  Polarization  data  of  a  GNF  electrode  (catalyst  loading  =  50  gg  cm  2)  in  0.5  M  H2S04  solution  in  the  absence  of  H202  (scan  rate  =  100  mV  s  ').  (b)  Determination  of  the 
apparent  total  double  layer  capacitance  by  plotting  double  layer  charging  current  (£  =  0.75  V  vs.  Ag/AgCl)  versus  scan  rate. 


1  to  0.5  V  vs.  Ag/AgCl  with  different  scan  rates  show  almost  a  flat 
current  plateau  suggesting  the  absence  of  faradic  reactions  in  this 
area  of  potentials.  Thus  the  potential  range  between  1  and  0.5  V  vs. 
Ag/AgCl  was  assumed  to  be  the  double  layer  region  and  the  current 
in  the  middle  of  this  region  was  employed  to  determine  the 
apparent  total  double  layer  capacitance.  Fig.  2(b)  presents  the 
measurement  of  the  total  double  layer  capacitance  of  a  GNF  elec¬ 
trode  at  E  =  0.75  V  vs.  Ag/AgCl. 

The  double  layer  capacitance  of  graphene  sheets  has  been  re¬ 
ported  between  2  and  21  pF  cnT2  [26,27],  By  assuming  first  a 
double  layer  capacitance  of  7  pF  citT2,  the  specific  active  surface 
area  of  GNF  was  165  ±  10  m2  g_1.  The  active  surface  area  deter¬ 
mined  by  voltammetry  is  less  than  50%  of  the  total  surface  area 
determined  by  the  BET  method  suggesting  that  most  of  the  catalyst 
surface  area  is  not  accessible  to  the  electrolyte.  If  one  assumes  a 
higher  double  layer  capacitance  reaching  30  pF  citT2,  the  specific 
active  surface  area  of  N-GNF32  is  found  to  be  36.7  ±  3  m2  g_1  (data 
not  shown  here). 

3.2.  HPRR  on  the  N-GNF32  electrode 

Fig.  3(a)  presents  the  steady  state  hydrodynamic  voltammo- 
grams  of  the  N-GNF32  electrode  in  1  mM  H2O2  +  0.1  M  Na2S04  at 
different  rotation  rates.  The  Koutecky— Levich  equation,  Eq.  (6),  can 
be  used  to  determine  the  apparent  number  of  electrons  n  involved 
in  the  HPRR  process,  and  the  apparent  heterogeneous  reaction  rate 
constant  k  of  H2O2  reduction  by  plotting  j-1  vs.  w-1^2  at  different 
electrode  potentials: 

1  _  _ 1  (6) 

j  nFkc  0.62nFDWaJ 


where  D  is  the  diffusion  coefficient  of  H202  (9  x  10~6  cm2  s_1) 
[22,28],  v  is  the  kinematic  viscosity  of  the  solution  (0.01  cm2  s_1),  c 
is  the  bulk  concentration  of  H2O2  (in  this  case, 
c  =  1  x  10~6  mol  cur3),  and  w  is  the  rotation  rate  of  the  electrode 
(rad  s_1).  Fig.  3(b)  shows  the  Koutecky— Levich  plots  for  the 
mixed  control  (mass  and  kinetic  control)  region  of  voltammograms. 
Straight  and  reasonably  parallel  lines  infer  that  the  reaction  is 
indeed  first  order  with  respect  to  the  H2O2  concentration. 

The  reaction  order  m  was  also  determined  as  1.16  ±  0.07  by 
plotting  log(j)  vs.  log(l  —  j/jd),  confirming  that  the  reaction  is  first 
order  with  respect  to  the  H2O2  concentration  (see  Fig.  SI  in  the 
Supporting  information). 

The  apparent  number  of  electrons  involved  in  HPRR  was 
determined  from  the  slope  of  the  Koutecky— Levich  plots  (Fig.  3(b)) 
at  different  electrode  potentials.  Fig.  4  shows  that  n  is  a  function  of 
the  electrode  potential  and  increases  while  the  electrode  potential 
shifts  toward  more  negative  potentials.  The  apparent  number  of 
electrons  less  than  2  suggests  that  the  decomposition  of  H2O2 
competes  with  the  reduction  of  H2O2  at  low  overpotentials.  While 
the  H2O2  decomposition  rate  is  independent  of  the  electrode  po¬ 
tential,  the  H2O2  reduction  rate  is  a  strong  function  of  overpotential 
and  increases  at  high  overpotentials.  Hence,  at  sufficiently  high 
electrode  potentials  ( E  <  -0.45  V)  the  H2O2  reduction  overcomes 
the  H2O2  decomposition  resulting  in  n  approaching  the  expected 
value  of  n  =  2  in  Eq.  (3). 

The  values  of  k  were  calculated  by  using  Eq.  (6)  and  the  in¬ 
tercepts  of  the  lines  in  Fig.  3(b)  at  different  potentials.  The  apparent 
heterogeneous  electron-transfer  rate  constants  strongly  depended 
on  the  electrode  potential  (Table  2),  suggesting  that  the  rate 
determining  step  is  an  electron  transfer  reaction  and  not  a  chemical 
reaction  [13,29,30], 


EPJ  vs.  Ag/AgCl 


Fig.  3.  (a)  Steady  state  hydrodynamic  voltammograms  of  the  N-GNF32  electrode  in 
1  mM  H202  +  0.1  M  Na2S04  at  different  rotation  rates:  100,  200,  400,  600,  900, 1600, 
and  2500  rpm  in  arrow  direction  (scan  rate  =  10  mV  s  ],  catalyst 
loading  =  250  pg  cnT2).  (b)  Koutecky-Levich  plots  for  reduction  of  1  mM  H202  at 
different  potentials:  -0.45  V  (□),  -0.5  V  (O),  -0.55  ( O ),  and  -0.6  V  (+).  Data  are 
obtained  from  Fig.  3(a). 


S 


-0.7  -0.6  -0.5  -0.4  -0.3  -0.2 


E/V  vs.  Ag/AgCl 


Fig.  4.  Apparent  number  of  electrons  for  the  reduction  of  1  mM  H202  on  the  N-GNF32 
electrode  at  different  potentials.  Data  are  obtained  from  the  slope  of  Koutecky-Levich 
plots  (Fig.  3(b)). 
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Table  2 

The  apparent  heterogeneous  electron-transfer  rate  constants  (k)  of  HPRR  on  the 
N-GNF32  electrode  at  different  electrode  potentials. 

£  (V)  -0.3  -0.35  -0.4  -0.45  -0.5  -0.55  -0.6 

k  (cm  s~’)  0.0023  0.0038  0.0056  0.0081  0.013  0.022  0.029 


3.2.1  Tafel  slope  analysis 

Tafel  equation,  Eq.  (7),  [31,32  was  used  to  construct  Tafel  plots 
at  different  rotation  rates. 


V  = 


2.3  RT 

acnaF 


(\ogj-  log  Jo) 


(7) 


where  17  is  the  overpotential  (V),  ac  is  the  cathodic  transfer  coeffi¬ 
cient,  na  is  the  total  number  of  electron(s)  involved  in  the  rate 
determining  step  and  all  previous  step(s),  and  jo  is  the  exchange 
current  density  based  on  the  geometric  surface  area  of  the  elec¬ 
trode  (A  cm-2).  Eq.  (7)  can  be  used  to  determine  the  ac  and  j0  from 
the  slope  and  the  intercept  of  the  line  between  r\  and  log(j). 

Fig.  5(a)  shows  that  for  high  enough  overpotentials  (higher  than 
35  mV)  two  Tafel  slopes  are  observed:  at  overpotential  more  pos¬ 
itive  than  -90  mV,  the  Tafel  slope  was  -127  mV,  while  it  was 
almost  doubled  at  overpotentials  more  negative  than  -90  mV. 
Doubled  Tafel  slope  phenomena  at  high  overpotentials  can  occur  in 
porous  electrodes  and  have  been  explained  in  detail  [33—35],  Tafel 
slope  analysis  was  repeated  for  different  rotation  rates  and  a  Tafel 
slope  of  125  ±  2  mV  was  obtained  suggesting  that  the  first  electron 
transfer  is  the  rate  determining  step  (na  =1).  Using  the  Tafel  slope 
at  different  rotation  rates,  na  =  1,  and  Eq.  (7),  the  cathodic  transfer 
coefficient  ac  can  be  evaluated  as  0.47  ±  0.01.  The  exchange  current 
density  based  on  the  geometric  surface  area  was  then  found  to  be 
(7.7  ±  0.3)  x  10-7  A  cm-2.  Using  the  active  specific  surface  area  of 
N-GNF32,  a  value  of  (8.3  ±  0.3)  x  1CT9  A  cm-2  was  determined  as 
the  exchange  current  density  based  on  the  active  surface  area  of  the 
electrode  (jew). 

The  mass  transfer  corrected  form  of  the  Tafel  equation  was  also 
used  to  determine  the  Tafel  slope  in  the  mixed  control  region  of 
hydrodynamic  voltammograms  [31,36,37],  Fig.  5(b)  shows  that  the 
Tafel  slope  at  different  rotation  rates  between  100  and  2500  rpm  is 
constant  (268  ±  2  mV)  and  independent  of  the  rotation  rate  sug¬ 
gesting  that  the  simultaneous  reduction  of  H2O2  and  O2  is  negli¬ 
gible  [13,38], 

3.2.2.  Stoichiometric  number  determination 

The  stoichiometric  number  (e)  is  defined  as  the  number  of  times 
the  rate  determining  step  is  occurring  in  order  to  complete  one 
cycle  of  the  overall  reaction.  The  stoichiometric  number  is  deter¬ 
mined  from  the  following  equation  [39]; 


Fig.  6.  Determination  of  the  slope  of  j  vs.  77  at  very  low  overpotentials,  (9j/<b?)?7— o-  Data 
are  obtained  from  Fig.  3(a)  for  w  =  100  rpm. 


nFj0 


(8) 


where  (9/794)ij-.o  is  the  slope  of  j  vs.  r\  when  jj  approaches  zero. 
Fig.  6  shows  the  relation  between  j  and  ij  at  very  low  overpotentials 
for  a  rotation  rate  of  100  rpm.  The  stoichiometric  number  was 
determined  to  be  1.14  by  using  the  slope  of  Fig.  6  and  Eq.  (8)  while 
jo  =  8.03  x  10-7  A  cm-2. 

The  stoichiometric  number  was  determined  for  different  rota¬ 
tion  rates  from  70  to  400  rpm  and  a  mean  value  of  e  =  1.17  ±  0.15 
was  obtained.  This  value  is  close  to  1  suggesting  that  the  rate 
determining  step  occurs  one  time  in  order  to  complete  one  cycle  of 
HPRR  on  this  catalyst.  The  stoichiometric  number  was  also  deter¬ 
mined  as  e  =  1  by  the  graphical  method  of  Allen  and  Hickling  (see 
Fig.  S2  in  the  Supporting  information)  indicating  consistency  in  the 
data. 

3.2.3.  FIPRR  mechanism  on  N-GNF32 

A  complete  mechanistic  analysis  of  HPRR  is  complex  and  its 
details  are  beyond  the  scope  of  the  present  work.  In  this  section,  the 
kinetic  parameters  determined  in  the  previous  sections  are  used  to 
propose  a  mechanism  for  HPRR. 

According  to  the  kinetic  parameters  determined  in  this  work 
and  different  mechanisms  proposed  for  HPRR  in  the  literature 
[13,40—42],  the  following  mechanism  is  suggested  for  HPRR  on 
N-GNF32  (NG  in  these  reactions  representing  the  nitrogen  doped 
graphene): 


NG+H202«NGH202 


(9) 


Fig.  5.  (a)  Tafel  plots  for  HPRR  on  the  N-GNF32  electrode  in  1  mM  H202  +  0.1  M  Na2S04  (cu  =  100  rpm).  Data  are  obtained  from  the  kinetic  control  region  of  Fig.  3(a).  (b)  Mass 
transfer  corrected  form  of  Tafel  plots  at  different  rotation  rates,  jd  is  determined  from  the  Levich  equation  by  using  n  =  2.  The  current  densities  are  obtained  from  mixed  control 
region  of  Fig.  3(a). 
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Fig.  7.  (a)  Hydrodynamic  voltammograms  of  reduction  of  1  mM  H2O2  in  0.1  M  Na2S04  solution  on  the  GNF  electrode  at  different  rotation  rates;  100,  200, 400,  and  1600  rpm  in  the 
arrow  direction  (scan  rate  =  10  mV  s-1,  catalyst  loading  =  250  jxg  cm-2),  (b)  Koutecky-Levich  plots  at  different  electrode  potentials;  0.1, 0.05, 0,  -0.05,  -0.1,  and  -0.15  V  vs.  Ag/AgCl 
from  the  bottom  to  the  top. 


NGH202  +  e-1^NG(0H)  +  0H-  (10) 

NG-(OH)  +  e-<-NG+OH“  (11) 

H202  +  2e~  ->20H~(The  overall  reaction)  (12) 

This  mechanism  satisfies  all  the  evidence  obtained  in  our  work: 
1)  the  first  electron  transfer  as  the  rate  determining  step,  2)  a  first 
order  reaction  with  respect  to  H202  concentration,  3)  a  stoichio¬ 
metric  number  equal  to  1,  and  4)  negligible  reduction  of  02. 
Moreover,  this  mechanism  is  consistent  with  the  mechanism  sug¬ 
gested  by  Wu  et  al.  [20]  for  H202  reduction  on  N-doped  graphene 
explained  in  the  introduction. 

At  low  overpotentials  both  H202  reduction  and  decomposition 
occur  on  the  electrode  surface  with  comparable  reaction  rates. 
Therefore,  reaction  (4)  competes  with  reaction  (10).  These  results 
support  the  equation  and  methodology  that  have  recently  been 
developed  in  our  group  for  analysis  of  HPRR  [43]. 

3.3.  HPRR  on  the  CNF  electrode 

In  this  section  the  H202  reduction  on  the  GNF  electrode  is 
studied  and  the  kinetic  parameters  are  determined.  In  order  to 
evaluate  the  catalytic  effect  of  nitrogen  on  the  graphene  struc¬ 
ture,  the  kinetic  parameters  of  HPRR  on  N-GNF32  and  GNF  are 
then  compared. 

Fig.  7(a)  presents  the  RDE  results  of  H202  reduction  on  the  GNF 
electrode  in  0.1  M  Na2SC>4  solution.  The  current  density  shows  a 
weak  dependence  to  the  rotation  rate  suggesting  that  H202  reaction 
rate  is  very  slow  on  the  GNF.  The  current  density  becomes 


Fig.  8.  The  relation  between  £  and  log(fc)  for  the  reduction  of  H202  on  GNF.  targe  Tafel 
slope  (-685  mV)  is  a  sign  for  an  electrochemical  reaction  with  a  slow  preceding 
chemical  reaction. 


independent  of  rotation  rate  at  £  <  -0.3  V.  These  results  suggest 
that  H202  is  relatively  stable  on  the  GNF. 

Fig.  7(b)  presents  the  Koutecky-Levich  plots  of  the  H202 
reduction  on  the  GNF  electrode.  The  apparent  number  of  electrons 
involved  in  the  HPRR  was  0.37  ±  0.07.  This  value  is  noticeably 
smaller  than  the  value  determined  for  N-GNF32  suggesting  that 
H202  decomposition  is  not  negligible  in  the  studied  potential 
window. 

The  apparent  heterogeneous  reaction  rate  constant  of  the  H202 
reduction  was  determined  from  the  intercepts  of  the  lines  in 
Fig.  7(b).  Generally,  the  reaction  rate  constant  of  an  electrochemical 
reaction  is  a  strong  function  of  overpotential,  and  a  linear  relation  is 
expected  between  r)  (or  £)  and  log  (k)  with  a  slope  equal  to  the  Tafel 
slope  [32], 

Fig.  8  presents  the  relation  between  £  and  log  ( k )  for  reduction  of 
1  mM  H202  on  the  GNF  electrode.  The  large  Tafel  slope  (-685  mV) 
observed  in  this  figure  can  be  a  sign  for  a  chemical  electrochemical 
reaction  mechanism  (CE  reaction  mechanism)  when  the  chemical 
reaction  is  the  rate  determining  step  [44].  This  result  suggests  that 
the  H202  decomposition  to  02  (Eq.  (4)),  followed  by  the  02  reduc¬ 
tion  (Eqs.  (1)  and  (2)),  is  not  negligible  on  GNF  in  the  mixed  control 
region  of  voltammograms.  The  CE  reaction  mechanism  can  also 
infer  to  the  slow  H202  adsorption  followed  by  the  H202  reduction. 

The  Tafel  slope  in  the  mixed  control  region  of  hydrodynamic 
voltammograms  was  determined  by  the  mass  transfer  corrected 
form  of  Tafel  equation.  Fig.  9  indicates  the  Tafel  slope  is  a  function  of 
rotation  rate,  confirming  that  both  H202  and  02  reduction  occur 
simultaneously  in  this  potential  region. 

The  Tafel  slope  and  cathodic  transfer  coefficient  in  the  kinetic 
control  region,  evaluated  using  Eq.  (7),  were  118  ±  2  mV  and 
0.5  ±  0.01,  respectively.  These  are  very  close  to  the  values 


Fig.  9.  The  mass  transfer  corrected  form  of  the  Tafel  slope  for  the  H202  reduction  on 
GNF.  Data  are  obtained  from  Fig.  7(a).  y  was  determined  from  the  Levich  equation  by 
using  n  =  0.45. 
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Table  3 

Kinetic  parameters  for  H202  reduction  on  GNF  and  N-GNF32  in  0.1  M  Na2S04  solution. 

n  Joa  (A  cnr2)  Tafel  slope  (low  overpotentials)  ac  e  F1202  decomposition  Simultaneous  reduction  of  02  and  H202 

GNF  0.37  ±  0.07  (1.52  ±  0.03)  x  10~9  118  ±2  0.5  ±  0.01  1  Yes  Yes 

N-GNF32  1.8  ±0.25  (8.3  ±  0.3)  x  IQ-9  125  ±2  0.47  ±  0.01  1  Yes  Negligible 


determined  as  125  ±  2  mV  and  0.47  ±  0.01  for  the  H2O2  reduction 
on  N-GNF32.  The  stoichiometric  number  was  also  determined  as 
1.12  ±  0.04.  These  results  highly  suggest  that  the  H202  reduction 
takes  place  by  the  same  mechanism  on  both  GNF  and  N-GNF. 

The  exchange  current  density  based  on  the  active  surface  area  of 
GNF.j'cm.  was  (1.52  ±0.03)  x  1CT9  A  cnT2  at  low  overpotentials.  This 
value  is  almost  6  times  smaller  than  the  determined  value  of  j'oa  for 
N-GNF32,  (8.3  ±  0.3)  x  1CT9  A  cnT2.  The  exchange  current  den¬ 
sities  based  on  the  active  surface  areas  indicate  that  higher  catalytic 
activity  of  N-GNF32  is  an  intrinsic  effect  of  nitrogen  atoms  on  the 
graphene  lattice,  irrespective  of  the  surface  area.  This  result  is 
consistent  with  the  result  reported  by  Wang  et  al.  [3]  for  02 
reduction  on  graphite  nitride.  They  showed  that  the  higher  catalytic 
activity  of  graphite  nitride  in  comparison  with  graphite  and 
graphite  oxide  is  created  by  nitrogen  atoms  which  facilitate  the 
electron  transfer  from  carbon  to  O2.  Moreover,  the  effect  of 
geometric  factors  such  as  surface  area  was  reported  to  be  negligible. 

A  summary  of  the  determined  kinetic  parameters  for  the 
reduction  of  H2O2  on  GNF  and  N-GNF32  is  presented  in  Table  3. 

3.4.  Comparison  between  HPRR  on  GNF,  N-GNF32  and  gold 

Fig.  10  presents  the  current  densities  produced  by  the  reduction 
of  1  mM  H2O2  in  0.1  M  Na2SC>4  solution  on  three  different 
electrodes:  GNF,  N-GNF32,  and  gold  at  400  rpm.  Comparison  be¬ 
tween  current  densities  produced  on  GNF  and  N-GNF32  indicates 
that  nitrogen  doping  increases  the  H2O2  reduction  rate  which  is 
consistent  with  the  determined  exchange  current  densities.  How¬ 
ever  the  produced  current  density  by  N-GNF32  is  still  lower  than 
gold  electrode.  Employing  N-doped  graphene  with  different  ni¬ 
trogen  composition  (higher  atomic  percentage  of  pyridinic  and/or 
graphitic  N  species)  and  using  other  dopants  such  as  iron  and  boron 
may  be  a  solution  to  increase  the  produce  current  density  in  order 
to  compete  with  gold. 

4.  Conclusion 

Although  both  catalysts,  GNF  and  N-GNF32,  employed  in  this 
work  present  high  BET  surface  areas,  determination  of  the  active 


EN  vs.  Ag/AgCl 


Fig.  10.  Hydrodynamic  voltammograms  for  the  reduction  of  1  mM  H202  in  0.1  M 
Na2S04  solution  on  three  different  electrodes  with  the  same  geometric  surface  area: 
GNF  (dashed  line),  N-GNF32  (dotted  line),  and  gold  (solid  line)  (w  =  400  rpm). 


surface  areas  show  that  more  than  50%  of  the  areas  are  not  acces¬ 
sible  for  the  electrochemical  reactions. 

The  results  indicate  that  N-doping  facilitates  the  H202  reduction 
on  graphene  while  it  does  not  change  the  mechanism  of  the  pro¬ 
cess.  Determination  of  the  Tafel  slope,  transfer  coefficient  and 
stoichiometric  number  at  low  overpotentials  indicate  that  H2O2  is 
reduced  by  the  same  mechanism  on  both  GNF  and  N-GNF32. 

The  H2O2  reduction  and  decomposition  occur  at  similar  reaction 
rates  on  N-GNF32  at  E  >  -0.45  V.  By  shifting  the  potential  toward 
more  negative  potentials  than  -0.45  V,  the  HPRR  becomes  domi¬ 
nant  and  the  number  of  electron  approaches  2.  Tafel  slope  analysis 
and  determination  of  the  apparent  heterogeneous  reaction  rate 
constant  in  the  mixed  control  region  of  voltammograms  reveal  that 
the  simultaneous  reduction  of  O2  and  H2O2  is  negligible  on  N- 
GNF32. 

Steady  state  hydrodynamic  voltammograms  demonstrate  that 
the  HPRR  occurs  at  limited  potentials  ( E  >  -0.3  V)  on  GNF.  At  po¬ 
tentials  more  negative  than  -0.3  V,  the  RDE  voltammograms  lose 
their  dependency  on  the  rotation  rate  suggesting  that  H2O2  is  stable 
in  this  potential  region. 

Although  N-doping  facilitates  the  reduction  reaction  of  H2O2, 
the  current  density  produced  on  N-GNF32  is  still  smaller  than  that 
on  a  gold  catalyst.  Further  studies  are  required  to  possibly  increase 
the  current  density  of  N-doped  graphene  catalyst  structures,  and 
understand  the  effect  of  different  N  types  involved  in  the  graphene 
lattice  or  from  other  dopants. 
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